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Abstract. We study the effect of Rashba spin-orbit interaction on the Josephson 
current through a double quantum dot in presence of Coulomb repulsion. In particular, 
' we describe the characteristic effects on the magnetic-field induced singlet-triplet 

£3 . transition in the molecular regime. Exploring the whole parameter space, we analyze 

the effects of the device asymmetry, the orientation of the applied magnetic field with 
respect to the spin-orbit interaction, and finite temperatures. We find that at finite 
temperatures the orthogonal component of the spin-orbit interaction exhibits a similar 
effect as the Coulomb interaction inducing the occurrence of a 7r-phase at particle-hole 
symmetry. This provides a new route to the experimental observability of the 7r-phase 
ON ' in multi-level quantum dots. 
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1. Introduction 

The Josephson effect is at the basis of Cooper pair transport through a weak link 
between two superconducting contacts. A current can flow due to the phase difference 
between the superconducting condensates; depending on the sign of the resulting 
current, one speaks of the 0- or the 7r-phase. The study of the Josephson effect 
is particularly insightful, if the junction itself has an internal structure, where the 
Josephson current is carried due to the formation of Andreev bound states. In quantum- 
dot structures a large tunability of the quantum-dot junction by electric gating or 
externally applied magnetic fields is provided. Recently quantum-dot superconductor 
hybrid structures (see Ref. [U [2] for a review) have been realized and the Josephson 
current as well as Andreev bound states have been studied in a systematic and controlled 
way |3liaiaEllZllSllSllinilIIllI2llI3lliaiiaiI5linilIE]. Whether the Josephson 
current through the system contributes in the 0- or the 7r-phase or whether it is 
completely suppressed, depends on the internal parameters of the quantum dot. In 
particular the Coulomb interaction has been shown to strongly affect the shape of 
the Josephson current [TUl SHI ETil EH E21 EH1 EH EEl EE], e.g. through the Kondo 
effect PIlSigHEllEBEEJEaffl Theoretical 
studies dealt with the single-level Anderson impurity coupled to BCS leads as a minimal 
model for the analysis of phase boundaries and the related transition in the Josephson 
current [22l [23l EH EU [32l [331 EH ESI [36l EZl SHI EHl SUl SU S2]- 

Various experimental realizations of Josephson multi-level quantum dots include 
carbon nanotubes, in addition to other materials such as InAs, both exhibiting strong 
spin-orbit (SO) effects [43J. While for the latter the relevance of the SO interaction is 
known [H] . for carbon nanotubes it has only lately been attributed to the curvature 
of the tube [I5j H6[ HZ]- SO coupling has an important impact on spin lifetimes in 
quantum dots J3H1H9}[50]. Several SO interaction effects on transport through quantum 
dots with normal leads have furthermore been studied, as e.g. the impact on weak 
localization and universal conductance fluctuations [51], the suppression of the Kondo 
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Figure 1. (Color online) Sketch of two serially coupled quantum dots with different 
energy levels ei/2 and onsite U, as well as interdot Coulomb interaction U'. The dots 
are coupled to each other by an interdot coupling with amplitude t, as well as by a 
spin-orbit interaction a. The external magnetic field B spans an angle 9 with respect 
to the SO direction. The double quantum dot is connected to superconducting leads, 
with gap Ae*' 8 , by a tunneling coupling t L / R . 
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ridges in combination with a Zeeman field [22], arid the occurrence of spin-polarized 
currents [53]. Observations of this type trigger the need to study the SO interaction on 
the Josephson current through interacting quantum dots. This is of particular relevance 
regarding the singlet character of the Cooper pairs carrying the Josephson current and 
the resulting sensitivity to interactions breaking the spin symmetry. Lately, different 
aspects of the SO interaction on the Josephson current and Andreev bound states have 
been considered in the regime of vanishing Coulomb interaction with an emphasis on 
the quasiparticle contributions to the current [251 EH EH1 EU [571 EEJ EHl EOj I6T] . 

In this manuscript, we focus on the effect of spin-orbit interaction on the Josephson 
current through quantum dots, allowing for a strong Coulomb interaction, where a 
magnetic field with arbitrary orientation with respect to the spin-orbit field can be 
applied. We thoroughly take into account the variety of parameters, governing a realistic 
quantum dot setup such as asymmetries in the coupling due to gating and contrast their 
properties with the characteristics of the SO interaction. We perform this study for zero 
as well as for finite temperatures. 

In order to assess the role of SO interactions in multi-level structures we here focus 
on a double quantum dot as a paradigm system, where the localized dots represent 
the different orbits. It has been shown in Refs. [211 122] , that the main features in the 
Josephson current through a single or a double dot setup can already be detected in the 
molecular regime of infinitely large superconducting gaps A. We therefore concentrate 
on the regime of A — > oo, which describes well the subgap features. In this regime single- 
particle transport is suppressed and the full spectrum can be accessed. Our detailed 
study can serve for the characterization of quantum-dot Josephson junctions and 
provides a reference for studies including quasiparticle transport as well [25l l26l l30 l 162]. 

The paper is organized as follows. In Sec. [21 we introduce the double-dot model and 
present the full effective Hamiltonian. In Sec. [3] we discuss the results for the Josephson 
current, starting from the symmetric double quantum dot, subsequently addressing 
the effects of the Coulomb interaction, SO interaction, finite asymmetries and finite 
temperatures. We finally extend this analysis to the more general case of a parallel 
quantum-dot configuration, where we focus on a special transport regime which cannot 
be accessed by the serial geometry. 

2. Serial quantum dot model 

2.1. Microscopic Hamiltonian 

We consider two serially coupled quantum dots with a single electronic orbital in each 
dot contributing to transport. All hopping amplitudes in the model are assumed to be 
real, except for the imaginary hopping amplitude ia taking into account the Rashba SO 
interaction. The model we use is sketched in Fig. CD The Hamiltonian of the isolated 
double dot is given by 

Hdd = Hq + H$ + Hi nt , (1) 
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containing the free part H , a term due to the SO coupling H SO} and the Coulomb 
interaction contributions, H int . The free part is 



b=i,2 



(2) 



with dj a being the creation operator of an electron on the dot site j = 1,2 with spin 
a =t, I- The on-site energies of the two quantum dots are 

B , & 

ei/v = -e + ± ^ > (3) 

which can be tuned by an external gate voltage and a magnetic field B with |B| = fi, 
lifting the spin degeneracy. (We here choose {iB = e = h~ = k B = 1.) The resulting 
Zeeman splitting is given by B with a = ± for spin up and down respectively; we neglect 
the effect of the magnetic field on the orbital motion. The Zeeman field B = Be z sets 
the spin quantization axis, which we here choose to point along the z-direction. The 
difference between the orbital on-site energies of the two dots is parametrized by the 
level detuning 6 = e^ a — e 2)CT . Hybridization of the single-dot states occurs due to the 
spin-independent interdot hopping with amplitude t. 

The Rashba SO interaction we are interested in is taken into account by an 
imaginary hopping amplitude of spin- dependent sign as the lattice realization resulting 
from spatial confinement in semiconductor structures [44^ 63J. Importantly, it also 
applies for the curvature-induced SO interaction which occurs in carbon nanotubes [4"5l 
l4*6l I4T] . The Rashba term of the Hamiltonian with a > reads 

H so = ia^2 [ d l,a( a z)a,a'd 1 vcos 9 + ^(^W'di^sin 9 + H.c. , (4) 

(7,(7' 

where 9 is the angle between the effective SO field a with |ot| = a and the Zeeman field 
B. The Pauli matrices a x and a z are related to the electron spin. The SO interaction 
can have a component parallel to B in ^-direction and a perpendicular one, which - 
without loss of generality - we choose to be parallel to the x-axis. For 9 = ±7r/2 the SO 
interaction and the 5-field are orthogonal, while they are parallel for 9 = 0. It can be 
shown that in the latter case the conventional hopping and the SO interaction can be 
combined to an effective hopping t = \/t 2 + 4a 2 . 

The local Coulomb repulsion between electrons occupying the double-dot is modeled 
as an on-site interaction U as well as a nearest-neighbor interaction U' 

H int = UJ2 (n jtt - (n u -^j+U' (m - 1) (n 2 - 1) (5) 

where the number of electrons on each dot is nj = J2a n j,o-i with = d- . Our 
model Hamiltonian is chosen such that e = corresponds to half filling. 

We consider the double-dot structure described above to be coupled to two 
superconducting leads. The left and right lead are modeled by the BCS Hamiltonian 

K = J' R = E e *4r<w - a J2 ( e ^ s 4t c U + H - c -) > ( 6 ) 

k,a k 
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with A and 4>l,r — ±0/2 being the BCS gap and phase, respectively. The coupling to 
the leads is modeled by the tunnel Hamiltonian 

H T = Y1 A* + *r4*A,* + H - c -] • ( 7 ) 

k,cr 

We furthermore assume the density of states p s in the leads in the normal state to be 
constant (this is particularly reasonable since we later assume large superconducting 
gaps). In this wide-band limit, the coupling strength to the leads is defined as 

r s = 7rt 2 sPs . (8) 

The total coupling is given by T = IY + Tr and we define the asymmetry in the coupling 
strengths to the two leads by (5 = (T L — r R )/T. 

2.2. Effective Hamiltonian for A — > oo 

We treat the problem by introducing an effective Hamiltonian in the limit where the 
superconducting gap is the largest energy scale. This allows to analytically determine 
the many-body eigenstates of the system as a starting point to exactly compute the 
Josephson current. 

We therefore set up the equations of motion for the double-dot retarded Green's 
functions Gf(t) = -i6(t){{%(t),*]{0)}), using the Hamiltonian introduced in the 

section above. Here, the so-called Nambu spinor of dot i is given by = (d^, d\,). The 
superconductor couples electrons and holes, while the particle number is not fixed, due 
to fluctuations in the superconducting condensate. Therefore it is useful to introduce 
the Nambu basis. This means that we count electrons of spin a and missing electrons of 
opposite spin a as Nambu particles. Importantly, the Nambu particle number is fixed 
in the presence of particle transfer due to Andreev reflection only. Due to the coupling 
to the superconducting leads also off-diagonal elements in this Green's function appear. 

In the limit A — > oo, the same set of equations of motion for G^f can be reproduced 
by an effective dot Hamiltonian, see the Appendix for details. For the double dot 
described in the previous section, this effective Hamiltonian is given by 

H cS = H + H so + H int + (T L e- t<t,/2 d lt d u + T R e llt)/2 d 2t d 2i + H.c.) . (9) 

Importantly, the proximity of the superconducting leads introduces a term which 
couples electrons and holes of opposite spin on each dot by Andreev reflection, but 
no term coupling particles on different dots is induced. We now individuate the 
relevant Nambu particle subspaces of the spectrum and their role for the Josephson 
current through the device. By introducing Nambu spinors, we can decompose the 16- 
dimensional Hilbert space of the proximized double quantum dot that underlies H c g 
into uncoupled sectors with different Nambu-particle numbers, {0, . . . , 4}. We note that 
the Coulomb interaction enters in a non-trivial way only the sector with Nambu-particle 
number 2. For the analysis of the ground state and the resulting Josephson current, it is 
helpful to relate the eigenstates of the effective Hamiltonian with a fixed Nambu-particle 
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number with the total spin s and the z-component s z . The Nambu-particle number is 
given by 2s z + 2. 

The characterization of the subspaces by the ^-component of the spin holds only 
in the absence of SO interaction or as long as the SO interaction is parallel to the 
magnetic field ck||S. For any finite orthogonal component a± with respect to the 
magnetic field, states with a difference \As z \ = 1 of the spin ^-component couple due 
to SO interaction-induced spin flips. The result is that only two independent subspaces 
are formed, characterized by an even, respectively an odd, Nambu-particle number. 
This coupling of different Nambu subspaces, will be shown to influence importantly the 
Josephson current. 

We are interested in the Josephson current through the device, which we can access 
through the full many-particle spectrum. Having calculated the set of eigenvalues Ei, 
the Josephson current at finite temperature can be obtained from the phase derivative 
of the free energy F [Ml [65] by 

J = 2<9 F, F = -Tln^V^ /T (10) 

i 

with F = — T In ^\ e~ Ei ^ T . This equation shows that in principle all eigenstates of 
the double quantum dot contribute to the total Josephson current with a contribution 
weighted by their energy. At zero temperature, T = 0, only the many-particle 
groundstate E contributes to the Josephson current. Therefore Eq. ( flQl) reduces to 

J = 2d < f ) E . (11) 

It depends on the parameters of the system, which subspace the many-particle ground 
state stems from. This can be tuned by a magnetic field B and the gate voltage entering 
the level positions €j. 

In the following section, we delineate the effect of the Coulomb interaction on the 
Josephson current in the double quantum dot introduced above, utilizing this as a basis 
to study in detail the effect of SO interaction and asymmetries. 

3. Results 

Based on the energy spectrum calculated from the effective Hamiltonian, Eq. (jHJ), 
we present results for the Josephson current in the following sections. We choose 
the coupling between the dots and between dots and leads to be of the same order 
t ~ r, except when otherwise indicated. In this regime, the eigenstates of the isolated 
double quantum dot are given by the molecular bonding and antibonding states, namely 
by the symmetric and antisymmetric superposition of the local quantum-dot states. 
Furthermore, we consider U = U' for simplicity, which is a reasonable choice for the 
double dot being in the molecular regime. 
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Figure 2. (Color online) (a) Josephson current for the symmetric (/3 = <5 = 0), 
interacting double dot with U — V as a function of the magnetic field B and the mean 
level position e in units of V. (b) Cuts of the density plot in (a) at fixed values of J5, 
as a function of e. The other parameters are t = T, k^T = 0, and 4> = tt/2. 

3.1. Effect of Coulomb interaction 

It has been shown before that for a symmetric double-dot setup without SO interaction 
and for sufficiently large Coulomb interaction, a singlet-triplet transition occurs as a 
function of an externally applied magnetic field resulting in a discontinuity in the 
Josephson current at a critical magnetic field [21]. We compare our results to this 
basic symmetric case and discuss in detail the impact of SO interaction, coupling and 
gate- voltage asymmetries and finite temperatures. 

We start our discussion by considering the effect of the Coulomb interaction on 
a fully symmetric dot, 5 = (3 = 0, in the absence of SO interaction (a = 0) and at 
zero temperature, k^T = 0. The Hamiltonian is symmetric with respect to B — > — B 
and we therefore restrict our analysis to B > 0. In Fig. [2} we show the result for the 
Josephson current through the double quantum dot, J(e, B), as a function of the mean 
level position, which can be tuned by applying a gate voltage, and as a function of 
an external magnetic field, at U = V . This generic scenario for the regime of large 
Coulomb interaction U > U c (with a critical interaction U c to be addressed later) has 
been discussed in Ref. [21 J. 

We show the Josephson current in units of a reference Josephson current at 
e = 0, which flows through a fully symmetric and non-interacting quantum dot, 
U = 5 = (3 = a = 0, at zero magnetic field and temperature, B = ksT = 0; for 
the plots shown in this paper it takes the value J (t = T, T, <\> = tt/2) = T/(2\ / 2 + V%), 
except when otherwise indicated. 

The Josephson current J(e, B) as displayed in Fig. [2] features three different phases, 
namely the 0-phase with J > (for < <fi < ir), the 7r-phase characterized by J < 
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(for < < 7r), and the triplet phase with J = 0. The respective ground state for 
the 0-phase is a non-degenerate singlet with {s = 0, s z = 0} (Nambu particle number 2). 
For the 7r-phase the ground state is associated with a free spin [TO, E2J EHl [SJ E7J EEJ 69J; 
at B = 0, this is given by the degenerate doublet with {s = 1/2, s z = ±1/2} (namely 
Nambu-particle number 1 or 3). The triplet state (Nambu-particle number or 4) has 
{s — 1, s z — ±1}, respectively. We furthermore note that J(e,B) is symmetric with 
respect to e — > — e. The symmetry in e is preserved as long as either 5 = or (3 = 
holds. 

In the vicinity of B ~ 0, depending on the gate voltage e, a positive Josephson 
current can flow through the double quantum dot in the 0-phase, associated with a non- 
degenerate spin singlet, see red regions in Fig. [2j However, due to the finite Coulomb 
interaction suppressing double occupation, a gate-voltage regime develops in which one 
of the molecular bonding and antibonding states is occupied with a single spin, leading 
to a reversal of the Josephson current. This latter is associated to a 7r-phase, where 
the double-dot state is in an almost twofold degenerate state with a free spin (exactly 
degenerate at B — 0). These two 0-7T-0 transitions as a function of e, occurring due to 
the single-level behavior, are observed as long as B < B c , where the critical magnetic 
field B c is defined by a degeneracy of a bonding and an antibonding state with opposite 
spin [22j[23l[2l[3lJ[32l[3^ this degeneracy corresponds 

to a degenerate singlet and triplet ground state. 

A finite Zeeman field B splits the spin degeneracy in the 7r-phase, leading at the 
same time to a broadening of the regions with negative Josephson current. The width 
of these regions depends on the strength of the coupling to the superconducting leads 
r and their phase difference <ft an d is of the order of U + B c for B > B c . The singlet 
remains unaffected by the magnetic field. For B < B c ~ t, the 7r-phase behavior is much 
more prominent than the singlet behavior close to zero gate voltage. Instead of a twofold 
0-7T-0 transition one would now rather speak of a single tt-0-tt transition as a function 
of the gate voltage. The resemblance of these features of opposite sign, at B ~ and 
B < B c , regarding the gate-voltage characteristics of the current J(e) as well as the 
current-phase relation J ((f)), can be perceived from panel (b) in Fig. [2] and is discussed 
in detail in Ref. |21j . Moreover, the line-shape of J(e) for B ~ B c is characterized by 
a discontinuity as a function of the magnetic field across the singlet-triplet transition 
at B c , as illustrated in Fig. [2j The suppression of J in the triplet spin configuration, 
where the occupation of the two levels with equal spin is favored, is due to the inhibited 
Cooper-pair tunneling with respect to the singlet one As a function of e, the width of 
the discontinuous singlet-triplet transition is of the order of the Coulomb interaction U 
for the symmetric case (as for the 7r-phase at B = 0). The occurrence of a singlet-triplet 
transition is obviously unique to a multi-orbital system. 

The discontinuity in J as a function of B at a critical magnetic field B c occurs only 

| This complete suppression of the Josephson current turns into a strongly reduced, but finite Josephson 
current in presence of quasi-particle transport occurring for a finite superconducting gap. However, 
these corrections turn out to be small as long as A ^> T. 
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Figure 3. (Color online) Critical Coulomb interaction U c as a function of the 
superconducting phase difference <f>. The other parameters are t — T, and kgT = 0. 
Inset: Critical SO coupling strength 2a c as a function of the superconducting phase 
difference <fr at U — 0. The other parameters are 9 = tt/2, t = T, and fc^T = 0. 



for sufficiently large Coulomb interaction ({7 > U c ), where the critical interaction U c 
depends on the superconducting phase difference, see the full green line in Fig. O We 
observe that for small values of 0, the singlet-triplet transition is observable already for 
rather weakly interacting quantum dots. For {7 = and <f> — > all three phases, 0, 7r 
and triplet, touch at e = and B = B c . In this case the discontinuous singlet-triplet 
transition as a function of B occurs hence only for e = 0. 

The main effect of the Coulomb interaction being smaller than the critical one, 
U < U c , can be observed in Fig. IH exemplified for the non-interacting, symmetric 
double dot. Note that the result remains qualitatively the same for any U < U c . The 
direct singlet-triplet transition, which for large U appears at B = B c , see Fig. [2j is 
inhibited by an interjacent 7r-phase, extending between the 0- and the triplet phase. 
The reason for this is that for weak Coulomb interaction and a magnetic field of the 
order of the single level spacing, a non-degenerate ground state with a magnetic moment 
associated to a spin 1/2 can exist at e ~ 0. The maximum strength of the Coulomb 
interaction allowing for the existence of such a spin- 1/2 ground state at e ~ is given 
by the critical interaction U c . 

We now focus on the S-field interval giving rise to the 7r-phase at e = 0, which 
prevents the direct transition between the singlet and the triplet phase. This interval is 
given by 



We will discuss the effect of SO interaction and of asymmetries in the following 
subsection. From Eq. (j!2p we see that for <fi — > the S-field interval, in which the 
7r-phase occurs at e = vanishes. This means that at (f> — )• a direct singlet-triplet 
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Figure 4. (Color online) Josephson current as a function of the magnetic field and the 
level position in units of T for U — (U < U c ). The other parameters are S = j3 = 0, 
t = T, k B T = 0, and = tt/2. 

transition as a function of -B occurs at e = already for vanishing Coulomb interaction. 
This finding is confirmed by the solid green line in the plot of Fig. [31 where the critical 
Coulomb interaction U c is zero at = 0. 

A further effect of the Coulomb interaction is that it tends to decrease the absolute 
value of the Josephson current. This can be seen from a comparison of Figs. [2] and HI It 
becomes even clearer in the quantitative comparison given by the solid green lines for 
the Josephson current as a function of the magnetic field at e = in the plots shown 
for U = and U = T, in Fig. El 

3.2. Spin-orbit interaction effects 

We now come to the main purpose of this paper, namely the study of SO effects. The SO 
interaction breaks spin-rotational symmetry by designating a certain (spin) direction. 
Spin is hence no longer a good quantum number, however at zero magnetic field a 
Kramers doublet remains as time-reversal symmetry is conserved. In the presence of an 
additional, externally applied finite Zeeman field, the direction of this external field with 
respect to the effective SO field is important. In the Hamiltonian responsible for the 
SO interaction in the double-dot setup, Eq. (jlj), the angle 6 between the SO interaction 
and the B field defines their respective orientation. A parallel configuration is similar 
to the situation discussed before. In this case, the amplitude of the SO interaction 
affects only the interdot tunneling, which takes a SO interaction-dependent effective 
value, t = y/t 2 + 4a 2 , and therefore the results of Sec. 13. II hold. 

In the following we focus on the effect of a finite orthogonal component a± and 
therefore concentrate on the case 9 = 7r/2 with a = a±. A finite orthogonal component 
of the SO interaction with respect to the Zeeman field introduces an anticrossing at finite 
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Figure 5. (Color online) Josephson current for the symmetric (5 = ^ = 0), interacting 
case, U — r, as a function of the magnetic field and the level position in units of T. The 
angle between the magnetic field and the SO interaction field is 9 = tt/2. Furthermore 
a = Q.5T, t = T, k B T = 0, and <j) = tt/2. 

B of the levels with a finite spin in either one of the dots or in both and consequently 
it leads to an avoided singlet-triplet transition, influencing the previously discussed 
discontinuity in the Josephson current as a function of B. For simplicity we consider the 
symmetric case with (3 = and 5 = and discuss the asymmetric case in the subsequent 
section. The effect of the SO interaction on the Josephson current is shown in Fig. |5l 
where we plot the Josephson current J(e, B) for a± = 0.5T and U > U c . The singlet- 
triplet transition as a function of B leading to a discontinuity in the Josephson current 
for a± = 0, as observed in Fig. [2j turns into a smooth crossover for a finite a± between 
the 0-phase and the triplet phase, where the Josephson current is strongly suppressed 
only in the limit of large magnetic fields. To show this effect quantitatively, we plot in 
Fig. E] (a) the Josephson current for e = as a function of B for various values of a±, 
clearly demonstrating the smooth crossover introduced by the perpendicular component 
of the SO interaction field. For comparison, in Fig. [6] (b) we show the corresponding 
non-interacting result, U = 0. For small SO interaction a direct singlet-triplet crossover 
is inhibited by an interjacent 7r-phase, see also Fig. HI However, for a± larger than some 
critical value a c , the 7r-phase is again not accessible at e ~ and the full (smoothened) 
crossover line between singlet and triplet phase is recovered. The Josephson current is 
hence positive for all B at e = 0. In this sense the SO interaction has a similar effect 
as the Coulomb interaction: it opens a window of positive Josephson current between 
the two 7r-phases, which appear due to transport when one of the different molecular 
states is occupied with a half-integer spin. This can also be understood by analyzing 
the relevant eigenstates of the system. The gap opened by the SO interaction in the 
single-particle energy spectrum of the isolated double dot leading to the anti-crossing 
is 2a± for U = 0. The critical SO coupling strength 2a c , is shown in the inset of 
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Fig. [3J The behavior is similar to the critical Coulomb interaction U c , however the 
slightly more linear functional dependence shows a modification of the spectrum due to 
SO coupling. At sufficiently large Coulomb interaction U > U c the width of the singlet- 
triplet transition as a function of e in absence of SO (in other words the minimal distance 
between the two 7r-phase contributions) is further enhanced by a finite perpendicular 
component of the SO field to U + 2a±. 

We point out that these effects inherent to the SO interaction depend on the angle 
between the SO interaction and the externally applied magnetic field. A rotation of the 
latter with respect to the SO field would therefore allow to address these effects explicitly. 
It is furthermore important to show that asymmetries between different dot states due 
to gating show features which are clearly distinguishable from the SO interaction. 
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Figure 6. (Color online) Josephson current for e = and different values of the SO 
coupling strength a as a function of the magnetic field in units of T. We show the 
symmetric case, 6 = /3 = 0, for a Coulomb interaction of (a) U = T and (b) U = 0. The 
angle between the magnetic field and the SO interaction field is 8 = n/2. Furthermore 
t = T, k B T = 0, and = tt/2. 
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Figure 7. (Color online) Josephson current for the asymmetric case, f3 = 0.7, 6 = T, 
as a function of the magnetic field and the mean level position in units of T for U = T 
(U > U c ) and = tt/2. The other parameters are t — T and ksT = 0. 



3.3. Asymmetry in coupling and detuning 

We consider, as a next step, the effects of asymmetries as generally present in 
experimental realizations of a double-dot setup on the features discussed before. The 
results are summarized in Figs. [3] and [7J 

The asymmetry due to different couplings to the leads, (3 ^ 0, and due to a level 
detuning, 5^0, result in changes in the 5-field interval, given for the symmetric case in 
Eq. ( 1T21) . and consequently in the critical Coulomb interaction, U c . Since the expression 
for SB in the asymmetric case is rather lengthy we do not present it here. However, 
we can extract the effect of various asymmetries from the dependence of the critical 
interaction U c on the superconducting phase difference 0, given in Fig. |3j A finite 
detuning 5^0 still allows for a direct singlet-triplet transition at zero interaction if 
4> — > and (3 = 0. Importantly, for a pure coupling asymmetry, (3 ^ and 8 = 0, a 
coexistence of the 0-, the ir- and the triplet-phase at zero Coulomb interaction U = is 
excluded and the required Coulomb interaction to allow for the singlet-triplet transition 
is strongly increased. This can be seen from the red dashed-dotted line in Fig. [3] In 
contrast, the figure shows that the overall effect of a finite detuning 5 is to reduce the 
critical Coulomb interaction U c . 

We now concentrate on the situation, where the Coulomb interaction is larger 
than the critical value at which the singlet-triplet transition occurs, and study the 
effect of asymmetries on the Josephson current. An overall effect of the two types of 
asymmetries, 5^0 and (3 ^ 0, is a reduction of the absolute value of the Josephson 
current. A fundamental difference in the discontinuity at the singlet-triplet transition 
arises furthermore when both asymmetry effects are present at the same time, namely 
a finite detuning and asymmetric dot-lead couplings, (3^0 and 5^0. In addition to a 
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suppression of the Josephson current, it exhibits a broken e — » — e symmetry, as shown 
in Fig. [7J For the strong asymmetry chosen here, (3 = 0.7 and 5 — t, a 7r-phase occurs 
at 5 = only for positive mean level energy e > 0. We can understand this in the 
following way: For (3 > 0, the left dot tends to be coupled more strongly than the right 
one; at the same time, for 5 > 0, the single-particle level of the right dot is lowered 
in energy and hence its overlap with the bonding state is increased. When the gate 
voltage is chosen such that e « y/S 2 + 1 2 /2, half of the energy difference between the 
single-particle bonding and antibonding level, the bonding state (mainly given by the 
right dot state due to the detuning) tends to be occupied with a free spin and Josephson 
coupling through the antibonding state is enhanced (due to (3 > 0). This leads to the 
7r-phase contribution at B — 0. However, when e ~ —\/ S 2 + 1 2 /2 the antibonding state 
tends to have a free spin, the transport channel through the bonding state is more 
weakly coupled ((3 > 0), and the 7r-phase only sets in at a finite 5-field. The effect, 
shown in Fig. [7J is reversed by changing either the sign of f3 or of S. 

This asymmetry goes along with yet a different effect, namely the tilting of the 
singlet-triplet transition line in correspondence of the discontinuity in J as a function 
of B, resulting in an e-dependent B c , differently to the symmetric situation. 

In the case where both the effect of SO interaction and asymmetries occur, the two 
effects superimpose. 



3.4- Finite temperature effects 

The analytic treatment of the effective Hamiltonian as outlined in Sec. I2.2[ allows 
to compute the exact Josephson current, J, not only for T = 0, but also for finite 
temperatures, see Eq. (flO|) . At finite temperatures, the transition between different 
phases is smeared out and the signal can be significantly reduced. Nevertheless, sign 
changes in the Josephson current have clearly been detected experimentally [HI EJ [10] . 
This reduction and smearing out of the signal due to finite temperatures is shown in 




Figure 8. (Color online) Josephson current at B = 0.2F as a function of e in units of 
the coupling strength V for different values of the temperature. The other parameters 
are t = T, = tt/2, U = T, f3 = 5 = and a = 0. 
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Figure 9. (Color online) Josephson current for e = as a function of the magnetic 
field in units of the coupling strength T for different values of the temperature. In (a) 
we set a — and examine the case of finite Coulomb interaction U = T. In (b) we 
show the non- interacting case U = for finite SO interaction a — 0.5T and 8 = ir/2. 
The other parameters are t — 2T, 6 = f3 = 0, and <f> = ir/2. The reference current is 
J o (£ = 2r,r,0 = vr/2). 



Fig. |8]for a = 0, where we show the Josephson current at B = 0.2T as a function of e 
for different temperatures (see Fig. [2] for the corresponding T = results). We clearly 
observe that the effect of the temperature is a general reduction of the current. Most 
importantly, effects from the contribution of higher energy states are expected to modify 
the shape of the Josephson current significantly. 

In the following we focus on the Josephson current as a function of the magnetic 
field at e = to monitor the impact of finite temperatures T > on the singlet-triplet 
transition (at a± = 0) and the smoothened crossover (at aj_ ^ 0). In Fig. [9J, we show two 
representative situations for the symmetric double-dot setup ((3 = 0,5 = 0), namely the 
one for finite Coulomb interaction, U > U c without SO interaction (a = 0) in the upper 
panel (a), and the one for finite orthogonal SO interaction, a± = 0.5T and a vanishing 
Coulomb interaction U = in the lower panel (b). At T = 0, for both parameter 
sets a direct singlet-triplet transition or respectively a crossover occurs, the 7r-phases 
of the Josephson current setting in only at finite values of e. The reason for this is 
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that the eigenstate of the system belonging to the subspace with one Nambu particle 
lies at higher energies for e = 0, and can therefore not contribute to the Josephson 
current if there is not a large enough temperature allowing for the thermal occupation 
of higher-lying states. 

We first discuss the effects at finite Coulomb interactions. In Fig. [9] (a) the 
discontinuity in the Josephson current J in correspondence of the singlet-triplet 
transition is gradually smeared out due to the increasing contribution of the higher 
energy states. Remarkably, in this regime the temperature stabilizes the 7r-phase, while 
in general the temperature leads to a reduction of the 7r-phase, in terms of a narrowing of 
the regions of negative current as well as in lowering the absolute value. This reduction 
of the 7r-phase at finite temperature is visualized in the Josephson current as a function 
of the mean level energy e, at B = 0.2r far from the singlet-triplet transition, in Fig. [HI 
where indeed the temperature smears out the features in the current leading finally to 
a complete suppression of the 7r-phase. However, in contrast to the effect of the SO 
coupling leading to a smooth crossover instead of a sharp singlet-triplet transition, the 
temperature smears out all transitions between different phases. 

A thermal stabilization effect of the 7r phase for a finite orthogonal SO interaction 
component with respect to the magnetic field a± = 0.5T is visible in Fig.[9](b), where the 
impact of the temperature on the Josephson current is reported for the non-interacting 
double quantum dot. The 7r-phase, which at e = and low temperatures is suppressed 
due to the SO interaction as in Fig.[5](b), is recovered for sufficiently high temperatures. 
The required temperature for this effect to occur depends on the distance between 
the energy of the lowest lying states of the even and the odd Nambu particle sector 
respectively, which in turn is directly related to the value of but also to the different 
tunnel-coupling strengths, t and T. In the situation shown in Fig. 0(b), the temperature 
recovering the 7r-phase is of the order of aj_. In this case, a sign change of an initially 
completely positive Josephson current J > at e = and <fi = ir/2 is obtained due to 
a temperature increase only. For a system in the absence of SO interaction, a± = 0, 
such an observation can only be expected in the presence of quasiparticle transport for 
a finite superconducting gap, where a small Josephson current can flow in the triplet 
phase. 

4. Parallel configuration 

Even though the physics of the serial double dot allows for a detailed insight in the 
physics of Josephson currents through quantum dots with SO interaction, also parallel 
quantum dots can be experimentally designed or can sometimes be a more appropriate 
model for a multilevel dot. We therefore consider in this section the situation of a double 
dot in the more general configuration containing also the parallel setup, as shown in 
Fig. EH with a finite coupling to both left and right superconducting lead for each dot, 




(13) 



l=L,Rj=l,2 a,k 
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We assume real and positive hopping amplitudes ty. The dot Hamiltonian and the 
lead Hamiltonian remain unchanged with respect to Eq. ([1]) and Eq. ([6]). In analogy 
to Sec. 12.21 (see also Appendix), we derive an effective Hamiltonian, starting from the 
equations of motion for the double-dot Green's functions in the limit of an infinite 
superconducting gap. Here, we find 

off = #dd (14) 

/2 N ' 



H. 



(r L1 e^ /2 + r me - 



(T L2 e irp/2 + r R2 e-^ /2 ) cLdt + H.c. 



n \ d ] rfi ■ 

For T L2 = T R1 = 0, the previous case of a serial double quantum dot is recovered. We 
observe that in contrast to the serial double dot, a parallel configuration allows also for 
crossed Andreev reflection [T7J [18j [701 [7H [72] [73] . The physics of the parallel geometry 
is very rich due to interference effects. We here refrain from giving a systematic account 
of the various transport regimes (see Ref. [74J for the classification of the different cases) 
and concentrate on one special situation which cannot be accessed by the serial geometry. 

In the following we address one of the cases in which destructive interference leads 
to the complete decoupling of one of the molecular levels. This is the case for symmetric 
couplings, T^i =Tl2 = = T_R2 = T/2, and vanishing detuning, 5 = [Til [751 [761 [77] . 
as can be seen performing a basis transformation of the effective (non-interacting) 
Hamiltonian to bonding and antibonding eigenstates of the isolated double dot 



\b,a) =_(|l )0 r> + |2,a» 
=^(|l,<7>-|2,a» 



(15) 



It turns out that the antibonding state is completely decoupled from the superconducting 
leads and the tunnel Hamiltonian in this basis reads 



iZt, 



. r ( e */a + e-W) 4 t dl + H.c. . 



(16) 



As a consequence, Cooper-pair transport is only possible through the bonding orbital. 

We start with the discussion of the non-interacting dot system, U = 0. Fig. [11] (a) 
shows the phase diagram of the many-particle ground state as a function of the gate 
voltage e and the magnetic field B. The three different possible types of ground states 




*L2 e 2,U t R2 



Figure 10. (Color online) Sketch of two parallel coupled quantum dots with 
parameters as in Fig. [T] 
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are, as discussed in the previous section, a singlet state with Nambu particle number 2 
(yellow), a triplet state with Nambu particle number (dark blue) and a spin 1/2 state 
with Nambu particle number 1 for spin —1/2 (green). The shape of the contributions of 
the spin 1/2 states reflect that the system is composed of one effective level coupled to the 
superconducting leads (giving rise to the rounded green region) and a second effective 
level which is uncoupled (leading to the sharp triangular feature). This behavior is 
confirmed by the results for the Josephson current, see Fig. [11] (b). For the Josephson 
current in this regime we find that indeed only the coupled bonding level contributes 
to the current. As a result we observe, that J is either positive, J > 0, or zero, J = 0. 
The reason for this is obviously that the spin of the antibonding level is not coupled 
to the bonding level in this highly symmetric system. Therefore the Josephson current 
does not reveal a 7r-phase as it was observed for the serial double dot in the single-level 
transport situation. The transport through the parallel dot is blocked, whenever the 
coupled bonding level is occupied with a spin 1/2 and hence Cooper pair transfer is not 
possible. 

As in the serial setup we use the situation of a finite Coulomb interaction as the 
reference situation to study SO interaction effects. The effect of Coulomb interaction 
is summarized in Fig. [12] (a-b). As for the non- interacting case, the decoupling of 
the antibonding state is reflected in a broken e — > — e symmetry. In presence of a finite 
Coulomb interaction, this broken symmetry in addition appears in a tilted singlet-triplet 
transition line of the ground state. The Josephson current shows that only a single level 
supports a Josephson current, while transport through the other level is blocked. Also 
the magnetic moment of the uncoupled level, when singly occupied, does not lead to a n- 
phase in the Josephson current through the other level. In contrast to the noninteracting 
case, the Josephson current does however change its amplitude when the ground state 
changes between a singlet and the doublet state, which can here be observed at negative 
e. The latter is due to a capacitive coupling between the two levels entering through 
the interaction term U. 

4-1. Spin-orbit interaction effects 

Also in the parallel setup the effects of the SO interaction can be clearly pointed out. 
We discuss the influence of the SO interaction on the Josephson current through the 
parallel double-dot setup, when the effective SO field is perpendicular to the externally 
applied magnetic field. The results are shown in Fig.[12](c-d). For simplicity we consider 
U — 0, the physical behavior is qualitatively not affected by a finite Coulomb interaction. 
Three different effects are visible due to the SO coupling: (i) it has a similar effect on the 
current as the Coulomb interaction in the sense that no state of spin ±1/2 is accessible 
in the vicinity of e ~ 0, allowing for a direct transition between a singlet and a triplet 
state, (ii) In contrast to the case without SO coupling, the discontinuous transition 
between positive and suppressed current is replaced by a smooth crossover as it was 
observed in the serial double-dot setup. This smooth transition again originates from 
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Figure 11. (Color online) Phase diagram of the groundstatc for a parallel geometry: 
singlet state (yellow), state with a free spin (green) and triplet state (dark blue) in 
(a), and Josephson current (b), for the symmetric, 8 = Pi = /?2 = 0, non-interacting 
case, U — 0, as a function of the magnetic field and the level position in units of T. 
Furthermore t — T, ksT — 0, and <j) = 7r /2- 

the coupling between the singlet and the triplet many-body ground state. This effect 
makes the regions in which the Josephson current tends to be blocked distinguishable. 
The reason is that in the subspace with one Nambu particle no current can flow through 
the bonding level when it is singly occupied. However, when the transporting bonding 
level tends to be blocked with one spin and at the same time also the antibonding level 
is, then SO coupling induces coupling to the current carrying singlet state. Therefore 
the region in which the current is fully blocked becomes much smaller (as expected for 
U = when comparing to the small 7r-phase in the noninteracting serial double dot 
studied before), (iii) We observe the absence of a 7r-phase originating from the fact that 
only a single level is effectively coupled to the superconducting leads being insensitive to 
the magnetic moment of the uncoupled level. The Josephson current however appears 
to be sensitive to the occupation of the (uncoupled) antibonding level through the SO 
coupling, which leads to discontinuities in the current, between the regions where the 
antibonding level tends to be occupied or empty. 
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Figure 12. (Color online) Density plot of the Josephson current for a parallel 
geometry, (a) for the interacting case, U = T and a = 0, and (c) for a finite SO 
coupling aj_ — 0.5T in the noninteracting case, {7 = 0. The current is shown as a 
function of the magnetic field and the level position in units of T. (b) and (d) cuts 
through the density plots in (a) and (c) of the Josephson current as a function of e and 
fixed values of B. The other parameters are the same as in Fig. [TTJ 



5. Conclusions 

We studied the effect of SO interaction a on the equilibrium supercurrent through a 
serially coupled double quantum dot as a minimal model to describe a multi-orbital 
system, attached to superconducting leads. The SO effects were illustrated in a study 
of the Josephson current and contrasted with effects due to Coulomb interaction and of 
coupling asymmetries. 

We find a pronounced dependence on the orientation of an applied magnetic 
field with respect to the spin-orbit interaction: for the SO interaction does 

not qualitatively affect the singlet-triplet transition. However, for a finite orthogonal 
component, a±, this transition is smoothened as the singlet-triplet transition is 
inhibited by an anticrossing of the corresponding levels. Extending the study to finite 
temperatures we find that a± has a similar effect as the Coulomb interaction in the 
stabilization of a 7r-phase arising around zero gate voltage. We finally discuss the 
Josephson current for the situation of an effective decoupling of double-dot states, 
which occurs for a symmetric parallel dot configuration, and point out the effect of 



Josephson current through interacting double quantum dots with spin-orbit coupling 21 

SO interaction in this more general setup. Our detailed analysis can serve as a basis 
for an extended study including quasi-particle transport due to a finite superconducting 
gap or through hybrid structures involving superconducting as well as normal leads, 
thereby allowing for an inclusion of non-equilibrium effects. 
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Appendix A. Derivation of the effective Hamiltonian of the serial double 
quantum dot 

The effective Hamiltonian can be derived by setting up an equation of motion for 
the full double-dot Green's function and showing the equivalence with the equation 
of motion of an effective Hamiltonian. The retarded Green's function, which we use for 
this purpose, is defined as G\jf(t) = — i9(t)({^fi(t), ^(0)}}, where the Nambu spinor 
tyj = {djp d^) has been introduced. The 6*-function assures that retarded processes are 
accounted for. Due to the vector form of the spinors, the Green's function is a 2 x 2 
matrix for electron and hole contributions. Off-diagonal elements appear due to the 
coupling of electrons and holes by the superconducting leads. We start by expressing 
the microscopic Hamiltonian of the full superconductor - double-dot - superconductor 
system, introduced in Sec. [2] in terms of Nambu spinors. At first we can omit the 
Coulomb interaction and the SO part for sake of clarity and show in the following that 
only terms in the Hamiltonian involving lead operators are affected by the coupling to 
the superconductor on the level of the Hamiltonian. Therefore the terms H mt and Hso 
are independent contributions which can simply be added to the effective Hamiltonian 
in the end of the procedure. Note however, that the results have to be modified if both 
dots couple to both leads as it is the case in the parallel setup. We start our evaluation 
from the microscopic Hamiltonian 

H = H + H lcad + H T (A.l) 

3=1,2 k,s 

+ + *tH TR % R + H.c.) 

k 

where = (c s jfcf , <4_ fe ,) is the spinor describing the lead operators and j equals 1 if 
j = 2 and vice versa. The Hamiltonian contributions in matrix form are explicitly given 
by 




(A.2) 
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We here performed a transformation, which absorbs the phase of the superconducting 
gap into the tunnel matrix elements, c s k a — > c^o-e^"/ 2 . 

The equation of motion is set up by considering the matrix equation i-^G^f(t) = 



5(t) - iO(t)({- [H,%(t)\ ,*i(0)|), see e.g. Ref. 
simplifies after a Fourier transform, with G 



jt ■- 



The resulting set of equations 
= / dte- iuJt Gji(t), yielding 



—H t co — H 2 — H R 



Gn G\2 
G21 G22 



1 . 



The expressions H s with s 
superconducting leads 



L, R introduced in Eq. 



(A.3) 

12) include the effect of the 



(A.4) 



They show that the A-dependence affects only those parts of the Hamiltonian which 
describe the coupling between the leads and the dots. Therefore we have to take the 
limit A —7- 00 only in H s . We here consider the wide-band limit, where the density of 
states in the leads in their normal conducting state is constant. We can then write the 
sum over the wave vector k as an integral and obtain 

Po |t,| 2 / co + e Ae^ 
*r - e 2 - A 2 1 Ae-^ s co - e 



de- 



VA 2 - co 2 



—CO 

-Ae- { 



Performing the limit A — > 00 results in 

-r s e^ s 
-r,e-^ s 



H. 



eff 



—CO 



(A.5) 



(A.6) 



The effective Hamiltonian in the limit of A — > 00 replaces the one given in Eq. flA.lj) 
and is given by 



H, 



eff 



E 



(A.7) 



Transforming Eq. (IA.7[) back into the electron basis and adding the Coulomb and SO 
interaction results in the effective Hamiltonian for the serial double quantum dot given 
in Eq. ([H]). 
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